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SECTION I - INTRODUCTION 


This report describes the physical and numerical formulation of a model for 
the horizontal solidification of a binary alloy. It can be applied, for 
example, to the case of steady-state solidification from two opposing walls 
in an ingot. The major purpose of the model is to calculate macrosegregation 
in a casting ingot which results from flow of I nterdendrf tic liquid during 
solidification. The flow, driven by solidification contractions and by gravity 
acting on density gradients in the interdendri tic liquid, is modeled as flow 
through a porous medium. Section 2 defines the symbols used throughout this 
report. Section 3 contains the physical formulation of the problem leading 
to a set of equations which can be used to obtain: (1) the pressure field, 

(2) the velocity field, (3) mass flew and (^) solute flow in the solid plus 
liquid zone during solidification. With these established, the model calculates 
macrosegregation after solidification is complete. Section 4 presents the 
numerical techniques used to obtain solution on a computational grid. Section 
5 contains results, evaluation of the results, and recommendations for future 
development of the model. Included are the macrosegregation and flow field 
predictions for tin-lead, aluminum-copper, and tin^bi smuth alloys as well as 
comparisons of some of the predictions with published predictions or with 
empirical data. 

The numerical model is the basis for algorithms used in a FORTRAN program which 
has been run on a Prime A00 computer system at Marshall Space Flight Center 
(KSFC) in Huntsville, Alabama. The program can be run in batch or interactive 
mode; the latter mode allows the user to interact with th'e model by changing 
process parameters or alloy, and by selecting graphical or tabular output for 
display on the terminal screen. The documentation for the FORTRAN code is 
Volume II of this report, and the operating manual is Volume III. Volumes I 
and III are self-contained in the sense that they describe the model and its 
operation without reference to the FORTRAN code; accordingly no programming 
or computer operation experience is needed to make use of the model. 
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SECTION 2 - NOMENCLATURE 


% m (Aj^jAy) 

8 

C" (wt. pet.) 

(wt. pet.) 
Cg (wt. pet.) 

Cg (wt. pet.) 
C s (wt. pet.) 
C Q (wt. pet.) 

g L . s s 

g (cm/s 2 ) 
g (cm/s 2 ) 
k 

K (cm 2 ) 

L (cm) 

M M 
n ssr SOR 

N i' N j 2 
p (dynes/cm ) 

0 (dynes/ cm ) 

z 

p Q (dynes/cm ) 
t (s) 

t £ , t L (s) 

T (°C) 
t e (°C) 
t l (°C) 

U (cm/s) 
t (cm/s) 

V x , Vy (cm/s) 
x (cm) 

Ax (cm) 
x E> x^ (cm) 
y (cm) 

Ay (cm) 


coefficient In the pressure equation, Eq. (3*4.3) 
coefficient In the pressure equation, Eq, (3-^-3) 
average composition defined by Eq. (3*2.4) 
composition of Interdendrl tic liquid 

local average composition of solid during solidification, 
defined by Eq. (3.2.5) 

composition of the solid at the sol Id- liquid Interface 

final local average composition, Eq, (3*5*0 

composition of liquid at the llquldus Isotherm, l.e., 
composition of the alloy 

volume fraction of eutectic liquid at the end of solidification 
volume fraction of liquid, solid during solidification 
gravitational acceleration 

gravitational acceleration In the -y direction 
equi 1 ibrium parti tion ratio 
permeabl 1 1 ty 
Ingot height 

maximum number of steady state or pressure Iterations 
number of mesh points in the x or y direction 
pressure 

modified pressure defined in Eq. (4.5.1) 

atmospheric pressure 

time 

time at passage of the eutectic front and liquidus isotherm, 
respectively 

temperature 

eutectic temperature 

liquidus temperature at composition C Q 

isotherm velocity 

velocity of the i nterdendrl t i c fluid 
x- and y~ t 'mponents of ^ 

horizontal distance from the chill face, see Figure 3-2 
grid spacing in the x- direction 

positions of the eutectic and liquidus isotherms, respectively 
distance from the bottom of the ingot 
grid spacing in y- direction 
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Y (cm 2 ) 

E SSI» E SOR 
P (g/(cm-s)) 

F (g/cnr*) 

P, (g/cnr ) 

L 3 

P s (g/cnr) 

P LE , P SE (g/cm 3 ) 

P LQ (g/cm 3 ) 


permeability coefficient 
convergence bounds; see Section 4.8 
vl scosi ty 

average density defined by Eq. (3.2.2) 
density of the i nterdendri ti c liquid 
density of the primary solid phase 

densities of the eutectic liquid and solid, respectively 
density of the liquid at the Hquidus Isotherm 


SECTION 3 - PHYSICAL DESCRIPTION OF THE MODEL 


3.1 PROBLEM DEFINITION 

Figure 3.1 shows an Ingot of a binary alloy undergoing horizontal solidification 
from chill faces located at opposite sides of an Insulated rectangular mold. 

The mold Is oriented so that the gravity force Is In the y-direction as shown, 
with no component normal to the x-y plane. Initially the contents are liquid 
with uniform composition C q . As heat Is extracted through the chill faces, 
the alloy solidifies bidirectionally from the chills and toward the centerline. 

It Is assumed that dendritic freezing takes place with a zone containing both 
solid and liquid phases (i.e., the S/L zone) which moves through the Ingot 
ahead of the eutectic Isotherms. Shrlnkage-and gravi ty -driven convection 
within these S/L zones causes non-unlformltles, known as macrosegregation, in 
the final local average composition after solidification, C s « The relationship 
between interdendri t 1c fluid flow and macrosegregation was first described and 
demonstrated in references 3~5 for flow due to solidification contraction only, 
and later for flow driven by gravity, in addition to solidification shrinkages 
in references 1, 2, and 6. The solidification model described below calculates 
the final local average composition after solidification by solving the hydro- 
dynamic equations for the flow of interdendri tic liquid within the limitations 
of several assumptions. 

The methods used here are similar to those of references 1 and 2 for macro- 
segregation resulting from steady-state solidification. The mold shown in 
Figure 3.1 is symmetric about the centerline, so it Is sufficient to treat 
only the left side of the ingot. Because there is no body force normal to the 
x-y plane, the flow in the S/L zone is two-dimensional. Nevertheless, much 
of the analytical development will be independent of the coordinate system and 
dimensionality of the S/L zone. As defined in reference 3, "unidirectional" 
solidification Implies unidirectional heat flow with planar isotherms, so that 
the S/L zone Is rectangular and the temperature gradient has no y component. 

In steady-state solidification the isotherms move through the ingot at a constant 
velocity, thus the geometry of the mushy zone and the temperature distribution 
within the S/L zone are constant in time. The macrosegregation resulting from 
steady-state solidification does not vary in the direction of heat flow; the 
final local average composition is a function only of vertical distance through 
the ingot. The 1 nterdendri ti c fluid flow is driven both by shrinkage due to 
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solidification contraction and liquid contraction during solidification, and 
by gravity acting on a fluid of variable density. The density gradient results 
from thermal gradients and more so from solute concentration gradients which 
must exist In the liquid within the 5/L zone. Flow resulting from the latter 
situation has been termed concentration-driven convection, and in references 1, 

2, 6 and 7 it has been shown to be a major contributing factor in the formation 
of channel-type segregates often called "freckles". A fundamental assumption 
of the solidification model Is no movement of the solid phase, so that the 
calculations cannot predict macrosegregation profiles for ingots In which 
freckles occur, although the conditions for and the location of the onset of 
freckle formation can be be predicted. The flow of interdendritic liquid is 
modeled as flow through a porous medium where the volume fraction available 
for flow is the local volume fraction of liquid during solidification, The 
movement of the S/L zone, the temperature field within the mushy zone, the liquid 
density as a function of composition and temperature, and the phase diagram 
(solid-liquid equilibria) for the binary alloy are assumed; the pressure field, 
velocity field, and macrosegregation are calculated. 

3.2 FLOW EQUATIONS 

Because the rate of isotherm movement is used as input to the solidification 
analysis and because no convection in the bulk liquid is assumed, the analysis 
can be restricted to the processes within the S/L zone and the solute flow at 
the boundaries of the S/L zone,, The coordinate system defined in the S/L zone 
is shown in Figure 3.2. The coordinates of any point (x, y) in the S/L zone 
are fixed relative to the ingot. The sides of the S/L zone at x^(t) and x^(t), 
the positions of the entectic isotherm and liquidus isotherms, respectively, 
move through the ingot with a constant velocity U. The analysis leading to a 
solidification model is based upon the concept of a volume element within the 
S/L zone which is small enough to be treated as a differential element, yet 
large enough so that the volume fraction solid within it is equal to the local 
average. The volume element concept is developed in reference 3. There is 
no solid phase movement into or out of the cell; there is no mass flux into 
or out of the element by diffusion; temperature and liquid density vary only 
differentially across the cell at any time. With these assumptions conversation 
of mass within the volume element is written 
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( 3 . 2 . 1 ) 


where 


§«- V.p L g L V, 
p - P s g s + P L g L . 


( 3 . 2 , 2 ) 


(See Section 2 for definitions of the symbols.) If solute enters or leaves the 
volume element only by liquid convection, If there Is no solute diffusion within 
the solid pha*:e, and If the composition of liquid within the volume element Is 
uniform,, then conservation of solute In the volu>«e element Is written 


^ (C p) - - 7.C L P L g L V, 
where C p - C s P s g $ + C L P L g L , 

and C s Is the local average solid composition at time t, defined by 

C s (x,y,t) - f p s C* (x,y,t) ~ (x,y,£)dt 

.."S J 

t, 


(3.2.3) 


(3.2.I|) 


(3.2.5) 


for t prior to passage of the eutectic front: t^<t<tg. 


Mass fluxes at the edges of the S/L zone provide boundary conditions for the 
conservation equations. The top and bottom of the S/L zone are In contact with 
the mold walls so that there Is no normal component of flow: 

V y » 0 at y * 0 and y * L for x E $xSx L> (3.2.6) 

As described in reference 1, flow at the eutectic isotherm must compensate for 
solidification shrinkage or expansion of the eutectic mixture: 


V ** 
x 


P SE~ P LE 
p LE ' 


U at x = x £ for 0<y<L. 


(3.2.7) 


Since It Is assumed that there is no convection in the bulk liquid, then 


V =0 at x = for 0<y<L. 


( 3 . 2 . 8 ) 


Additional conditions at the liquidus isotherm are 
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(3.2.9) 


C L » C Q and ■ I at x ■ x L for 0£y£L, 

Darcy's Law for flow through a porous medium Is used to describe the flow of 
the Interdendrl ti c liquid In the S/L zone; thus 


$ - • (?p + p l 3), ( 3.2.10) 

L 

where p Is the viscosity of the liquid and K Is the permeability given here, 
as In references 1 and 2, by 

K-Y9l (3.2.11) 

with y a constant. This form of the permeability function Is isotropic; J.e., 
there Is no directional preference In the resistance to flow. 


In addition to the assumptions of steady-state, unidirectional solidification, 
which lead to planar Isotherms moving with a constant velocity, a constant 
horizontal temperature gradient Is assumed. Then the temperature field in the 
mushy zone is given by 


t.-t, 

T(x,y,t) - T £ + J^(x~x E ) 


( 3 . 2 . 12 ) 


for x^x^x^ and 0<y<L, as shown in Figure 3.4. The steady-state assumption 
now gives the Isotherm velocity 


U 


3T / 31 „ 3T / / T l" T E, 
9t ' 9x 9t ' 


(3.2.13) 


Equations (3.2.12) and (3.2.13) show that the temperature field and isotherm 
movement are determined upon specification of the cooling rate, 9T/9t, and 
the mushy zone width, (x^-x^) . The liquid in the volume element is assumed 
to be in equilibrium with the solid at the solid-liquid interface so that the 
phase diagram relates liquid composition to temperature in the volume element. 
As shown in Figure 3 • 3 » the liquidus of the phase diagram is approximated by 
a straight line, and the composition of the interface solid is given by 

C* = kC L , (3.2.14) 
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where k Is the equilibrium partition ratio. C L and 9C L /9t can thus be calculated 
as functions of temperature by 


dC. 

C L (T) * C E + dT 


(T-T e ) and 3t 


BC. dC 


L 9 T 
d'T 9t 


(3.2.15) 


where Cg, T £ and dC^/dT are taken from the alloy phase diagram. The resultant 
linear variation of C L through the S/t zone is shown in Figure 3. i*. Liquid 
density is a function of temperature and composition, but the equilibrium 
assumption allows density to be written as a function of composition only, 
Re?>, ences 1 and 2 have shown that the density of the interdendrl tl c liquid 
can be approximated by a linear function of composition and that the density 
of the primary solid phase is approximately constant In the solidification 
temperature range. Thus liquid density and 9p L /8t can be calculated by 


dp. 9p 

p L< C L> " P LE + dC^ (C L~ C E> *" d 3t 


dp L 3 C l 
dC[" STt ' 


( 3 . 2 . 16 ) 


P s can be considered constant, and the densities of the solid and liquid phases 
are completely specified by known values of p s , p $E , dp L /dC L , and p LE> The 
variation of p L across the S/L zone is shown in Figure 3.1< where p L increases 
during solidification. Some alloys have a decrease of liquid density during 
solidification and they can be treated, as well, providing the correct sign 
for dp^/dC^ is used. 


An equation used to relate g g and g^ is provided by the assumption of no pore 
formation: 

9g s 9s, 

9 s + 9 l * ' , ° r at a — • (3.2.17) 


Equations (3.2.1) through (3.2.17) and the associated assumptions completely 
describe the solidification model. However, the solution for the velocity 
field and local fraction liquid can be simplified if the information above is 
used to derive two more equations which replace the continuity equations in 
the solution process. These equations are derived in Sections 3.3 and 3.k. 
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3.3 LOCAL SOLUTE REDISTRIBUTION EQUATION 


The mass and solute conservation equations can be combined to get an equation, 
derived originally in reference 3, that gives the local fraction liquid In 
terms of the local velocity and the local liquid composition. By expanding 
the right side of Equation (3*2.3) and substituting Equations (3.2.1) and 
(3.2.2) , we arri ve at 


3g, 


-y.C L P L g L tf « cJL (p L g L ) - C L P S “- - P L gJ*VC L . 


(3.3.1) 


When Equation (3*2.5) is differentiated, there results 


3t ^ 


S r jl r * 
g s at °s g$ u s at 


1 39c * 

i“ at” (c s- c s ) 


(3.3.2) 


Expanding the left side of Equation (3*2.3) and substituting Equations (3.2.14), 
(3.2.17) and (3.3.2) leads to 


9_ 39 

fit 


ac 


“1-1 a 

(c p)= " kc L p s at" + p l 9 l at - + c l 3 T^ p l 9 l^ * (3.3*3) 


Finally, Equations (3-3.1) and (3*3.3) are substituted into Equation (3.2.3), 
giving 


39, 


3C. 


1_ Zk - (J-\ JL rllk . \ 

g L at " i-k p $ c L ( at v V * 


(3.3.4) 


Since the liquid composition and temperature are directly related (see Equation 
3.2.15), then 


ac. ac. ac. r. T 

VC L= ar VT and ar FT fp 


(3.5.5) 


and Equation (3.3.4) becomes 


9 l 3t 


(_i_) (l + 2 ^ 1 , 

' 1-k' p_ U fiT 1 c. at • 
* at L 


(3.3.6) 


Equation (3*3.6) is called the local solute redistribution equation. 
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3*4 PRESSURE EQUATION 


Numerical .solution of the sol fdl ft cation model equations can be facilitated by 
replacing the mass conservation equation with an elliptic equation for pressure. 
Using Equations (3.2.2) and (3.2.17) and the assumption of constant solid density 
to expand the left side of Equation (3.2.1) gives: 


9.P 

at 


9g, Bp. 

^ p l' p s^ at - + 9 l at“ * 


(3.4.1) 


The pressure equation results from substituting Equations (3.2.10) and (3.4.1) 
I nto Equation (3.2.1): 

fkp, *1 3g. 3p, 

V- [— (Vp+p L g)J- (P L -P S ) ~ 9 L ^~- 0 , (3.4.2) 

Equation (3-4.2) can be written In the form: 



v 2 p 

+ ft*Vp 

+ B 

* 0 



(3.4.3) 

wi th 

ft. 

_,H_ y 
Kp L 

Kp 

(__L) 

M 

f 



(3.4.4) 



— H-f v 

kp L L 

Kp, 

.V 

Sg. 

8p l1 

(3.4.5) 

and 

B = 

•( L 
y 

p L g) - 

^ p l' p s^ W 

‘ " 9 l atj 

As discussed in 

Section 

4.5, 

several 

1 numerical 

techniques 

are avai lable for 

solution 

of equations w 

i th the form 

of Equation 

(3.4.3). 



Conditions (3.2.6) through (3.2.8) can be rewritten to provide boundary conditions 
on pressure at all sides of the S/L zone. Substituting Equation (3.2.6) for 
in Darcy's Law, Equation (3.2.10), gives 


8p 

97 = ■ p L g 


and substituting Equation 


3p = ^SE 
3x “ K p 


at y = 0 and y = L for x^<x<x^, 
(3.2.7) for V in Equation (3.2.10) 

X 

~ P LE 

U at x = Xp for 0£y<L. 

LE 


gives 


(3.4.6) 


(3.4.7) 
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If there Is no convection In the bulk liquid, then 

p ■ p + p, 9 (U-y) at x ■ x, for OSySL. (3.^.8) 

0 L o L 

3.5 MACR05EGREGATI0N 


After solidification Is complete 
average composition Is given by 


C s (x,y) 



at a (x,y) In a casting, 

Cg(x,y,t)dg s (x,y,t) + P SE g E C £ 

P s (l _9 E ) + Pse 9 E 


the final local 


(3.5.1) 


where the llquldus Isotherm passes the point (x,y) at time t^ and the eutectic 
front passes at t c . The integral term accounts for microsegregation within the 

U 

dendritic solid and the second term in the numerator accounts for the solute in 
the solidified eutectic mixture. 

For steady-state solidification C,. is a function only of y, and the integral 
in Equation 3.5.1 is equal to the same integral taken between the limits x ( to 
x £ at a fixed time t. Macrosegregation as a function of height in the ingot 
is given by c‘ s (y)-C o . 

3.6 SUMMARY 

The equations which constitute a complete description of the solidification 
model implemented in the computer program are as follows: 


Local Solute Redistribution Equation 

1 39 L / 1 n P L T i , ^-VT ~| L_ 

i[3t =- W P S |I J C L 3t 


( 3 . 6 . 1 ) 


with g L = 1 at x = x L for 0<y^L. 


( 3 . 6 . 2 ) 
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I 


Pressure Equation 

V 2 p + t-V p + B « 0 


for 
and 
wl th 


Kp, 


(Zk), 

n Mr. ¥ ' tl ' » 


Kp 


11 r Kp L + 9 9| 

k^L V,( 1T p L g) " ^ P l" P S^ 3T" 9 l 

—• - -p L g at y * 0 and y * L for x^SxSx^ 

II at x * x,, for 0<y<L, 
P LE fc 


3p 


3p Pg L P SE~ P LE 


and 


3x K 

p “ p o + P L 9 ^ L " y ^ at x “ \ for OSySL 


Porous Flow Model 


* 


u - - ctr (Vp + p l 5) 


K » Y9l 


Temperature Field 


V T E , v 

T “ t e + (x - x e ) ,or x e SxSx l 


'L E 


u = - 2E/3I 

3t / 3x 


Liquid Composi tion 

C L “ C E + dT~ (T -V ’ 9t ~ dT at 


dc. ac L dc L 3T 


Liquid Density 


_ n dp L , v 8p L dp L dC L 3T 

P LE dC, VU L V * 3t “ dC. dT 3t 


Final Composition 


c s (y) 



kc L dg s + P SE g E ^ E 

^(l-g E ) + P SE g E 


(3.6.15) 


The following data are required as input to the solidification model. 

dC L 

Phase Diagram : T^, C^, k 

Alloy ; C q 

Solidification Process : L, (x^-x^) , 9T/3t, g 

Density; p L£ , dp L /dC L , p $ , p gE 
Vis cos i ty : y 

Permeability Coefficient : y 

Upon specification of these values, the description of a steady-state casting 
of a binary alloy is complete, and Equations (3-6.1) through (3.6.15) can be 
solved to yield the final local average composition as a function of vertical 
distance through the ingot. The following section describes the numerical 
analysis which is the basis of a computer program which performs such a 
calculation. 


SECTION k - NUMERICAL METHODS 


The equations presented In Section 3 require reformatting In a discrete form 
for solution on a digital computer. Numerical techniques applied to the 
continuous equations yield finite difference equations and procedures for 
calculation of their solution. This process Involves selection of finite 
difference methods that are appropriate to the type of partial differential 
equation to be solved as well as systematic verification of the convergence of 
the discrete solution to the solution of the continuous equations. Although 
the last fifteen years have seen tremendous advances In the discipline of 
numerical analysis applied to fluid flow equations, there Is as yet no standard 
approach to the process of solving a particular set of equations. Techniques 
that provide a valid solution to one problem can produce an Invalid but 
physically plausible result when applied to a slightly different problem. For 
this reason the verification phase of the process Is crucial: a "solution" 

should never be accepted merely because it has a qualitative physical explanation. 
Procedures for sufficient validation of a solution to a flow problem are even 
less routine than the selection of the solution technique; they require experience 
In recognizing the behavior of an invalid result as well as an understanding of 
the foundations of tlWi numerical approximation process. 

Beyond the necessity of finding numerical methods which produce a sufficiently 
accurate solution, the constraints of interactive computing on a time-shared 
mini-computer require the use of memory-efficient, rapidly-convergent techniques. 
Numerical methods applied in the context of batch computing on a large high-speed 
mainframe are frequently chosen for their ease of application and because of their 
familiarity; in an interactive mini-computer application such methods can cause 
excessively long response times. Thus the use of extremely small grid spacings 
and double precision arithmetic, standard "fixes" for inadequate numerical methods, 
are to be avoided in min i -computer applications in favor of numerical techniques 
that are more specifically tailored to the problem at hand. 

The techniques described in this section have been implemented and validated in 
an interactive program on the Prime ^00 at MSFC. Under heavy job-load conditions 
the response time for the primary calculation phase is less than five minutes. 


i»-l 


4.1 COMPUTATIONAL GRID 


The numerical solution Is calculated on a rectangular grid covering the S/L 
zone as shown In Figure 4.1, The value of a function at a grid point Is 
considered representative of the behavior of the function In the surrounding 
cell so that, for convergent numerical methods, the solution of the finite 
difference equations is close to the solution of the partial differential 
equations when the mesh spacing is sufficiently small. In the following 
discussion subscript notation indicates location in the mesh: for any 

function defined in the S/L zone, 

f jj - f (*j »Yj ) 

where Xj * x^ + (i-l)Ax for i * 1, ...» Nj , 

and Yj * (j-l)Ay for j = 1, ..., Nj . 

Ax and Ay are given by 

■ ispr and ^ = yt • 


4.2 OVERVIEW OF SOLUTION PROCESS 


The outline of the solution process is the same as the one used in reference 2. 
Details of each step may vary from those employed in reference 2 because of the 
difference in coordinate systems as well as the use of alternate numerical 
techniques. The solution steps are explained in detail in later subsections 
and the program modules containing the relevant code are indicated. 


1 . 
2 . 

3- 

4. 


Read the program input as defined in subsection 3-6. 


9C. 


3p, 


Calculate X, U, C^, -g^-, p L and -g^— at each point in the S/L zone by 
evaluating equations (3.6.11) through (3-6.14). 


3g, 


using the analyti 


Calculate an initial estimated distribution of g L and g-^- 
solution for the zero gravity case, equations (4. 3- 1 ) , (4. 3- 2) , and (3-6.1); 
evaluate K according to equation (3*6.10). 

Calculate an initial estimated pressure solution from equations (4.5.12) 
and (4.3-2), 

Evaluate the pressure equation coefficients, pT and B as shown in subsection 
4.4. 


5 - 







6 . 

7. 

8 . 


9. 


10 . 


Solve the pressure equation using the relaxation technique discussed In 
subsection *1.5. 

Use Darcy's Law as shown In subsection *1.6 to evaluate the velocity of the 
Interdendrl tic flow. 

Use the velocity calculated at step 7 In the local solute redistribution 

equation, as shown In subsection *».7, to calculate new distributions of 

9g. 

g L and ; re-evaluate K by applying equation (3*6.10) at each point In 
the mesh. 

Check convergence of the nonlinear system by applying the criteria described 
In subsection *1.8. If It has not converged and If freckling has not been 
detected, repeat steps 5 through 9. 

Calculate the final local average composition. 


*1.3 INITIAL ESTIMATED SOLUTION 

Starting the iterative calculation of the steady-state solution requires an 

9g. 

estimate of g, and •*— - to be used In generating the pressure equation coefficients 

L. (7 V 

for the first pass through steps 5 through 9 . Although the Initial estimated 
solution has no effect on the converged solution, convergence to the steady- 
state g L , if one exists, can be speeded by starting the iteration with a good 
estimate of g^. Hence the analytic solution for the case with no gravity force 
is used to estimate g^. The only driving force In this case Is shrinkage 
contraction so that the y-component of velocity Is zero. The expressions for 
g L , derived in reference 1, is: 


C|> <ir-'-i)c E 

K se 


a 


g L - 


C o a + <r- -'>**“ 

se 


(*K3.1) 


wi th 

9g 


a = - 


1-k * 


3t 


— can be. evaluated directly from equation (3.6.1) with Vy « 0 and 


V - - 
x 


gi (P s - p L ) + g E (p SE -p s > 


p L g L 


*,_*, 


u. 


'(* 1 . 3 . 2 ) 


9 * ■» 


The Initial estimated solution Is calculated In subroutine EST. 
h.k EVALUATION OF THE PRESSURE EQUATION COEFFICIENTS 

For general distributions of there Is no analytic form, so the gradient and 
divergence terms In and B must be evaluated numerically. Comparison of centered 
finite difference formulas with evaluations of the analytic expressions available 
In the ze:o gravity case shown that differencing the logarithmic form of the 
spatial derivative terms Is an order of magnitude more accurate than differencing 
the forms shown In equations (3.6. l|) and (3.6.5). Using the logarithmic form 
provides a more accurate solution for a given mesh size, or It permits a coarser 
mesh In achieving a given accuracy, thus reducing computation time. The 
logarithmic forms of A and B subject to g * (o,g) are: 


3g L 

9t~ 


V [ ln( -y L_) ] ’ 

and B - gp L fy [ ,n( ir">] “ [ (p L" p s } 

9p L 1 

+ g L ir J • 

Centered finite differences applied to these forms give: 

_ r Kp. Kp. *1 

( Vij [ ln ( ~ ) i+i ,j “ ,n ( T” ) i-i,j J /(2Ax) 

for l<i<Nj and lSjSNj, 

. r Kp, Kp, 1 

( V i j “ [ ln ,n ( "ir , i,J-lJ /(2Ay) 


{hJ ♦.!) 


(4.4.2) 


(4.4.3) 


(4.4.4) 


for l<i<N. ( //nd ]<j <Nj , 


and 


8_ 

ay 


r 2 i r 2 ~i 

Kp* ' * Kp . 1 Kp /■ 

1 n (“ )J i j = ) i , j + 1 ~ 1 n { ~ ] f , j - 1_ 


/ (2Ay) (4.4.5) 


for 1£i£N. and Kj<N.. 

' J 


At the edges of the mesh centered differences cannot be used; the x-dcrivative 
terms are approximated by first order difference formulas of the form 
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(4.4.6) 


(fl) « t (f^ j-fjj)/Ax for 1*1 , N j and l<j<Nj, 


and the y-derlvatlve terms are approximated by second order difference formulas 
of the form 


(|~) - + I3 f |j-'» f f j+, +f ( j+ 2 )/(2Ay) for J-l f Nj and 1<1<N ( . (4.4.7) 


Subroutine PSETUP evaluates A and B. 


4.5 SOLUTION OF THE PRESSURE EQUATION 


Because atmospheric pressure, p Q , Is very large compared to the pressure 
variations In the S/L zone, the accuracy of a numerical solution for total 
pressure would be considerably less than that of a solution for p-p Q . Similarly, 
the large hydrostatic pressure near the bottom of the S/L zone causes an 
unacceptable loss of significance in the calculation of (p-p Q ) during the 
the evaluation of V . These inaccuracies, Inherent in direct solution for p, 

A 

can be avoided by solving a modified pressure equation for 

P * P-P Q - P Lo g(L-y) (4.5.0 


where p (Q g(L-y) is the bulk fluid hydrostatic pressure. f5 is the solution of 

V Z p + rf'Vp + B = 0 (4.5.2) 

***V 

where A Is unchanged and 

B = B-p Lo gA y . (4.5.3) 

The boundary conditions for equation (4.5.2) are 


3P 

9y 

3p 

3y 


(p Lo “Pl^ 9 at V~° and y =L for x L <x<x E , 


W l P sr p LE 

K Ple 


U at x = x £ 


for 0<y£L, 


(4.5.4) 
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(4.5.5) 


(4*5,6) 


{> » 0 at x * x^ for O-y-L 

Standard solution procedures exist for linear elliptic equations such as 
equation (4,5.2). Iterative techniques have been preferred for a number of 
years, although fast direct (non-iterative) methods have been developed recently 
for problems with a large number of mesh points (references 8 and 9), Direct 
methods can be cumbersome to apply, especially If there are derivative boundary 
conditions or first order terms In the elliptic equation. Iterative techniques 
tend to converge slowly, but they are applicable to the most general form of the 
elliptic problem. The solidification model requires only a relatively small 
number of mesh points, It has derivative conditions at three boundaries, and It 
has a first order pressure term. Consequently, the Iterative technique known 
as successive overrelaxation (SOR), or the accelerated Gauss-Setdel scheme, was 
applied to the pressure equation, and the strategy developed In reference 10 
was used to speed convergence by dynamic estimation of the optimum acceleration 
parameter. 


The recursive formula for the (n+l)st Iterate of fS according to the method of 
successive overrelaxation (reference 8, p. 468) Is 


p^j* 1 * (l-w)5 


U 


(4.5.7) 


, , . n ~ n+1 n n+1 o 1 

* - P M,J) * + + 6 B ijJ 


for l£i£Nj and lSjSNj 


where 


V 


9 * T x - , 

x 2(Ax fc + Ay z ) 

a 2 

Ax 


(“xhj ■ 


6.. “ 9 9 , 

y 2 (Ax + Ay ) 

2 2 

$2 _ Ax Ay 


“ 2 Ay ^Vij’ 


2 (Ax 2 + Ay 2 ) 


and 0 ) is the relaxation parameter described below. The restriction i<N. is the 
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boundary condition at x - x^: P N j does not change from Its Initial value, 0. 

At the other boundaries the offse; terms are evaluated according to the formulas: 



a n+1 

S p 2 j - 2Ax(||) j j 

at i ■ 1, 

(4.5.8) 


a n+1 
HJ-1 

1 P i2 " 2Ay( f^li 

at j - 1, 

(4.5.9) 

and 

« n 
^ij+1 

+2Ay # ) IN j 

at J * N. , 

(4.5.10) 


where the partial derivatives are given by the boundary conditions, equations 
(^•5*4) through (4.5.6). This method of applying the derivative boundary conditions 
is second-order accurate in Ax and Ay as is the centered differencing used at 
Interior points. A more popular derivative boundary technique which calculates 
the boundary values directly as, for example, 


,n+l 

} i! 


ft n-H 

P , 2 




cannot be applied to the present problem because It would introduce an inconsistency 
at the corners (x^, 0) and (x^, L) where a normal derivative condition meets o 
pressure condition. Although the latter technique is easier to apply and has a 
heuristical ly pleasing form, Its use would lead to significant numerical error 
in V near the corners. 

r> 


The strategy developed by B. A. Carre® in reference 10 is used to estimate the 
optimum value of the relaxation parameter w. Use of the optimum w was observed 
to reduce the number of iterations in the SOR technique by a factor of ten to 
fifteen over use of u=>l , the Gauss-Se’^el scheme. The optimal w is estimated 
dynamically during solution for p as follows: 


1. Choose an initial w in the range l<w<u> 0 where is the (unknown) optimal 
value of a). 

2. Perform N iterations on equation (4.5.7). 

on set 
N, N. 

I 

l i = l 


After the iteration set 


r (n, = 1 t k ( , n) - e I }"- ,) | 

j=l i=l * * 

for n = N- 1 and n == N, 
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3. The optimums is approximately 


4. 


5,. 


6 . 


'V/ 

CO * 

o 



A+ ut_- 1 

, ~1 
AO) J 


i 



Set to * to Q - ^(2-to Q ) 

Repeat steps 2 through I» until W Q stops changing: 
Continue Iterating on equation (I/,5.7) with to=to Q 


Ato 

o 

2 -to 


<.05. 


until p converges. 


The values of w and N in steps 1 and 2 are taken as to=’. 375 and N«12. 


The convergence test on p checks the largest relative change In p over each 
iteration until 


max 



ft (n-l) 





+ 




(A. 5. 11) 


where e_ B is controlled by the progress of the steady-state solution as 
described In section * 1 . 8 . If condition (I».5.H) is not obtained within 
iterations, an error message is displayed and the calculation proceeds to step 

7- 


An initial p distribution is needed to start the SOR iterations. Although the 
initial p does not affect the final pressure solution, a substantial savings in 
computing time can be achieved by starting with a good initial estimate. During 
the first iteration on the steady-state solution the analytic zero-gravity solution 
is used to start the pressure iterations; for subsequent steady-state iterations 
the final p distribution from the previous steady-state iteration is used. In 
the zero-gravity case V^=0, so j5 can be calculated by solving Darcy's Law, 

equation (It. 6 . 1 ), for J|- and integrating it across the S/L zone: 

x V9, , 

P (x,y) = 0 + J (--^V x )dx 5. 12) 

x L 

V for the zero gravity case is given by equation (It. 3 . 2 ). 


All procedures described in this subsection are coded in subroutine PSOLVE except 
the initial estimated p which is calculated in subroutine EST. 


4,6 VELOCITY CALCULATION 

After the pressure distribution is known, the velocity can be evaluated directly 
from Darcy's Law, equation (3.6.9), written in terms of |5: 

* - ™tvp + (p L -P Lo )g). (^6.1) 

Ha L 

At interior points the gradient term is replaced by centered finite differences: 
(ff>ij - for KKN. and lSjSNj, 

and (|i) fj - for 1<I<N. and l<J<N j . 

First order one-sided differences of the form (4.4.6) are used at x^ while the 
boundary conditions provide values for the remaining pressure gradient terms* 


Velocity is calculated in subroutine VLCTY: K is evaluated in subroutine PERM. 


4.7 SOLUTION OF THE LOCAL SOLUTE REDISTRIBUTION EQUATION 


After the velocity field is known g^ can be calculated by solving the local 
solute redistribution equation. In the steady-state case this can be simplified 
to an integration along lines of constant y. The steady-state hypothesis, equation 
(3*5.2), applied to equation (3-6.1) gives 


J_^L 

g L sx 



Then g^ is approximated by 

(g. ) . . = (g, ) . , , . 
a L 1 j 3 L i+l ,j 


<C L>ij 


-q 




+1,J 


(4.7.1) 


(4.7.2) 


with (g L } N. 


' ,J 


q = i (q ; j 


i + l ,j 


), 
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V, - * 

s 


and q j j defined as the quantity in brackets in equation 1) . After g L is 

known, Sg^/St is calculated by a direct evaluation of equation 

g^ and 3g^/3t are calculated In subroutine LFRAC, 

4.8 ITERATION TO THE STEADY-STATE SOLUTION 

At the end of each iteration on the steady-state solution, steps 5 through 9 
in Section 4.2, g^ is tested for convergence. The convergence criterion is a 
comparison of the maximum relative change with an input threshold: 


0 (m) (m-i)| ,i/ (m) (m-l)\L F 

g L 9 U | /i(g L + s L )J< E ssi 

> j 


max 
i 


(4.8.1) 


where g^is g^ at the end of the m^ iteration. If condition (*1.8. 1) is not 
obtained by the end of M ssf iterations an error message, "steady-state solution 
did not converge" is put on the screen and the iteration process is terminated. 
E, 


'SSI 


and Mgjjj are user-specified inputs as described in Volume II! 


The convergence threshold for the pressure iteration, © S0 R’ is contr °^ ec * by 
the progress of the steady-state iteration. Initially the pressure convergence 
threshold is 10 E SQR , when the left side of inequality (A. 8.1) drops below 100 
Eg S | » e sOR ' S rec * ucec * to 5 ^sOR' an< ^ w ^ en t ^ )e s '^ e °f (4.8.1) is below 10 

E or , , , e cnD is equal to E SQR . This strategy speeds up the calculation by avoiding 


'SSI’ SOR 

unnecessarily tight convergence of the pressure equation while g^ is relatively 
far from its steady-state value. 


Iteration to the steady-state solution is controlled by subroutine SSICON. 

*4.9 CALCULATION OF MACROSEGREGAT i ON 

The integral in equation (3.6.15) can be transformed to an integral with respect 
to the logarithm of g^: 

*E X L rl 

f kC L dg g = J' kC L dg L = J kC L g L d(ln g L ) 

X L X E X E 

* 4-11 


(4.9.1) 


Accurate numerical evaluation of the second Integral In (4.9.1) would require 
an extremely fine mesh for alloy systems such as Sn-Bl whose large k causes g^ 
to be steep near the llquldus Isotherm* The last Integral In (4.9.1) Is amendab! 
to accurate approximation on a relatively small mesh. It Is approximated as 

X L N I 

/ kC L g L d(ln g L > - E ^kipC^), + p" 9^-ln 9 L| _J> t'l.9-2) 

* 

the remaining terms In (3*6*15) are evaluated directly, with g^ equal to g^ 
at i=*l. 


C Is evaluated in subroutine MACSEG. 
s 

^.10 FRECKLE CONDITION 


The existence of freckling in a given case can be detected by the program. 


I I ic ii cur\ i lily Cui iui uuu i o i uun i i cUte i l my . 


|I = |I + tf*VT>0 (4. 10.1) 

When condition (4.10.1) is detected at any point in the mesh, the calculation 
is stopped. 


The freckle condition .is tested in subroutine FRECKL. 


SECTION 5 - EVALUATION 


This section presents and describes the results of calculations using the model 
for flow of interdendri ti c liquid and macrosegregatlon described in Sections 3 
and I*. Calculations are for three alloys (Al-^.5% Cu, Sn-15$ Pb, and Sn-3% Bi) 
which show the effects of varying (l) the height of the solid-liquid zone, (2) 
the gravity force, and (3) the cooling rate during solidification. Added 
discussion Is given to the formation of "freckles" and to the differences in 
macrosegregation behavior among the three alloys. 

Ingot Height 

Figure 3.2 illustrates the solid-liquid (S/L) zone in a unidirectional ly 
solidified horizontal ingot. The width of the S/L zone is - X^, and the 
height is L. Figure 5.1 shows the effect of ingot height on macrosegregation. 
As ingot height increases, the major effect is that segregation is reduced, 
particularly in the central portion of the ingot. Thus the taller the ingot, 
the greater is the portion of the ingot near C D (the average composition of 
the aiioy). in this regard, Figure 5.1 compares well with Figure 10 In 
reference 1. However, there are differences at the top and bottom of the 
ingot which are attributed to differences in the choice of finite difference 
equations used to formulate the boundary conditions at the top and bottom of 
the ingot and in the method of solving for the velocity of the i nterdendr i t i c 
liquid. In this work, since velocity is sensitive to the volume fraction of 

liquid, then a scheme is used whereby an iteration calculation is done to solv'a 
for velocity. Algorithms used herein give results with improved accuracy of 
those given in Reference 1; details are given in Section A of this report. 

Gravity Force 


Figure 5.2 shows the effect of the gravity force on macrosegregation in Al-^.5% 
CU alloy. With zero gravity, there is no macrosegregation because flow of 
i nterdendri t i c liquid is solely due to solidification shrinkage. As such, all 
flow is normal to isotherms as shown in Figure 5*3a. The effect of gravity 
is to cause the flow to be downward within much of the S/L zone, and there is 
some reversal of flow in the lower right areas of Figures 5.3b and c. Notice, 
also, that the magnitude of the vectors increases with increasing gravity force. 
Figure 5.^ gives the solute flux corresponding to Figure 5.3 in terms of g of 
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Figure 5.1* Effect of ingot height on macrosegregation. 





effect of gravity force on Macrosegregation in 
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t.M* 


(b) 


tin* 


» MMUI 


(c) 


Figure 5.3. Effect of gravity on the velocity of i' nterdendri ti c liquid in 
A1 - k.5% Cu: (a) 0 g; (b) i g; (c) 2 g. Conditions are the 

same as in Figure 5.2. 









'Z’S sJnBtj u| sb siues sip sjb suoitjipuoo 
•63 (o) ‘6( (q) ‘60 (ej :Aonv no %$-ij - iv u i *p jnbi i 
o i q j jpuspjsqu i sip ui xny sqn^os sip uo A 4 |abj 6 rjos^-jo 


-fj'S sjn6 1 j 


i mi • M-mtt i 


I UMI -Wl I M-tHlI i — 

•*14 MW 




t I I 1 t 1 « • I • > I > • 

f t t I t 1 t t t t t t f • 

r t t t r t t » t t t t t t 

/ r r i t t t t t t t r t t 

/ / l ! t t t f t t t t t t 

/ ff/rtrifttvit t 
/ / / /./ t t t t t t t f r 
f/frrrtttffrri 
//////•ftfrrrft 

SSSSS/ttftffe i 

/ / / / * 4 * * * 

4 * 4 4 • 


«• 

0 



a 




• 

« 

• 

• 


9 

/ 

9 

f 

9 

9 

« 

• 

• 

• 

* 

*—• 

9 

9 

9 

/ 

9 

9 

9 

9 

9 

9 

• 

— 

9 

9 

9 

/ 

/ 

t 

9 

9 

9 

9 

9 


00 

9 

rf 

/ 

/ 

/ 

9 

9 

9 

9 

9 

m00 

*9 

9 

9 

4* 

/ 

/ 

9 

9 

9 

9 

9 


00 

/ 

s 

/ 

/ 

/ 

/ 

if 

9 

9 

* 


09 

✓ 

S 

s 

/ 

/ 

/ 

/ 

9 

9 

9 

— 

00 

/ 

s 

s 

/ 

s 

/ 

/ 

9 

9 

9 


09 

09 

* 

s 

/ 

9 

9 

9 

9 

9 

9 

— 

00 

s 

s 

s 

✓ 

4» 

9 

9 

9 

9 

9 


00 

0~ 

s 

s 

4* 

✓ 

9 

9 

9 

9 

9 


90 

s> 

9» 

** 

4* 

✓ 

9 

9 

9 

9 

9 


09 

9* 

** 

** 


** 

0 * 

9 

9 

9 

9 

— 

90 

s9 

09 


0” 


0 * 

9 

9 

9 

9 

— 

000 


00* 


0+ 


00 

00 

00 

*0 

*0 


__ 




_ 














2 

Cu/s-cm , These fields are similar to those of Figure 5*2 and are given here 
to point out that this option for output Is available. Figures 5.3 and 5.4 
also show why segregation Is negative at the top and positive at the bottom 
In Figures 5.1 and 5.2 when a gravity force acts on the liquid. Since the 
density of the liquid Increases during solidification, the solute-rich liquid 
flows away from the upper portion towards the lower portion of the Ingot 
(Figures 5*3 and 5.4) resulting In negative segregation at the top and positive 
segregation at the bottom. For alloys In which liquid density decreases during 
solidification, flow would be from bottom to top (opposite to Figures 5.3 and 
5.4), and segregation profiles would be reversed from those shown In Figure 5.2. 

Cooling Rate 

For a given gravity force, the most Important process variable, which influences 
macrosegregation, is cooling rate. This Is shown In Figure 5.5 In which 
macrosegregation curves for Ig at three different cooling rates have been 
drawn. One of the curves is that given in Figure 5*2 with a cooling rate of 
-0.333°C/S. When the cooling rate Is doubled (-0. 667°C/S) , the segregation is 
substantially reduced. However, when the cooling rate is reduced to -0.100°C/S, 
segregation is increased significantly. 

The effect of cooling rate on segregation can be explained by comparing the 
magnitudes of the velocity of i nterdendrl tic liquid and the Isotherm velocities. 

With slow solidification, isotherms move slowly and there is more flow through 
a unit volume, within the mushy zone, during the time for complete solidification; 
thus macrosegregation is increased. Figure 5.6 shows this effect for the two 
extreme cooling rates of Figure 5.5. With the slower cooling rate (Figure 5.6a), 
the magnitude of the velocity of the interdendri tic liquid is less than that at 
the greater cooling rate (Figure 5.6b), but relative to the velocity of the 
isotherms, there is substantially more flow at the lesser cooling rate. Consequently, 
macrosegregation increases with decreasing cooling rate as shown in Figure 5.5. 

Formation of Freckles 


When flow of interdendritic liquid is extensive and macrosegregation is severe, 
there is a reversal of that flow. This can be seen in the lower right areas of 
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Effect of cooling rate on the velocity of i nterdendri t ic 
liquid in A1 - *».5% Cu: (a) -0.100°C/s, isotherm velocity 

is 5.05 x 10"3 cm/s; (b) -0.667°C/s, isotherm velocity is 
3.37 x 10“ 2 cm/s. Conditions are the same as in Figure 5.5 




Figures 5.3c, 5.*»c and 5.6a. Ther* can be cases, which have been documented by 
experiment, In which the flow of the InterdendrI tic liquid In the direction of 
the Isotherms, exceeds the Isotherm velocity. When this occurs, there is local 
remelting of the dendritic solid. The remelted region Is the source of an 
unstable channel which persists within the solid-liquid zone and results In a 
severe localized segregate (I.e., "freckle 11 ) which Is rich In eutectic constituent 
after solidification Is complete. The criterion for this phenomenon Is 

($.VT/e)-i (5.1) 

where V Is the velocity of the InterdendrI ti c liquid, VT is the temperature 
gradient and e Is the cooling rate. An alternative expression for equation 
(5.1) Is 

(n.\?)/(n*tl) £l (5.2) 

In which n Is a unit vector normal to the isotherm and U is the velocity of 
the Isotherm. Equation (5.2) essentially states that If the velocity of the 
liquid, In the direction of the moving Isotherm, exceeds the velocity of the 
isotherm, then local remelting occurs. Equation (5.1) or equation (5.2) as an 
equality can be used to predict the onset of tno formation of a freckle. However, 
if the inequality applies, then the model developed herein does not accurately 
predict macrosegregation. 

In Figure 5.6a, the maximum value of the L.H.S of equation (5*2) is approximately 
0.3; therefore, no freckles are predicted. For this particular Ingot, the onset 
of freckling occurs at a cooling rate between -0.033°C/S and -0.100°C/S. 

Pressure Field 

As discussed in Sections 3 and * 1 , calculation of the pressure in the solid-liquid 
zone is an intermediate step to compute the velocity of the interdendri tic liquid. 
The major purpose of the analysis of macrosegregation is to calculate macrosegre- 
gation, as in Figures 5-1, 5-2 and 5.5; complementary information sought is often 
the velocity field of the interdendri tic liquid as in Figures 5.3 and 5-6. There 
are circumstances, however, in which the pressure field would be of interest; one 
example would be to predict evolution of dissolved gas during solidification. 
Figures 5-7 (a) and (b) shows the pressure field corresponding to the velocity 
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fields of Figures 5.6 (a) and (b), respectively. In Figures 5.7 (a) and (b) , 

P Q Is the ambient pressure at the top of the liquid pool (approximately 1 atm 
or 10^ dynes/cm 2 ) and P Js the pressure In the solid-liquid zone. Except at the 
bottom of the ingot solidified with a cooling rate of -0.100°C/s, there Is a 
decrease of pressure through the sol id- 11 quid zone. This pressure drop, of 
course, is attributed to the friction associated with flow of liquid through 
the dendritic network. As cooling rate Increases, or permeability decreases, 
the pressure- drop increases. For example, at y/L * 0.5> the pressure drop 
across the mushy zone for a cooling rate of -0.100°C/s is about 600 dynes/cm 2 
(Figure 5*7a); at an increased cooling rate of -0.667°C/s, the pressure drop 
across the solid-liquid zone at y/L = 0.5 is about 2800 dynes/cm . 

Other A) loys 

Figures 5.8 and 5-9 show results for flow of interdendri tic liquid and macro- 
segregation, respectively, in Sn-I5% Pb alloy in an ingot of similar size as 
described above for Al-^.5% Cu alloy. in the Sn-Pb alloy, the volume fraction 
of liquid at the eutectic isotherm is approximately 3 -A times that in the Al- 
Cu alloy (0.35 vs. 0.09), so that the average volume fraction of liquid in the 
solid-liquid zone is also significantly greater in the Sn-Pb alloy than in the 
Al-Cu alloy. As a result, flow is very extensive in the Sn-Pb alloy (Figure 
5.8) and much more so than shown in Figure 5.6. Hence, macrosr "Vjgation is 
also extensive as demonstrated by Figure 5-9 where there is a spread of 9% Pb 
from top to bottom of the ingot. Since the dendritic network of this ingot is 
so "open", then the pressure drop through the solid-liquid Is almost non- 
existent; this is seen in Figure 5.10. 

Finally, $n-3% Bi alloy is an example of an alloy with a relatively minor amount 
of liquid (approximately 1-2%) at the eutectic isotherm. Also, it dif ers from 
Al-*».5% Cu alloy and Sn-15% Pb alloy in that there is an expansion on solidifi- 
cation near the end of freezing. Thus in Figure 5.11, the magnitude of the 
velocity vectors is about one order of magnitude less than in Figure 5.8, and 
the effect of expansion (rather than contraction) is apparent in the left side 
of the sol id- liquid zone. Macrosegregation is only about 0.5% Bi from top to 
bottom in this ingot (Figure 5.12). The pressure field, Figure 5.13, is interesting 
in that the pressure "drop" is reversed corresponding to the direction of flow 
for much of the sol id- liquid zone of Figure 5.11. 
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Figure 5-11. The velocity of interdendriti c liquid 
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APPENDIX A 

ALLOY PROPERTIES AND SOLIDIFICATION PARAMETERS 



WEIGHT PERCENT COPPER 









DENSITY (g/ cm 5 ) TEMPERATURE (°C) 











P SE (density of eutectic solid) « 8.60 g/cm^ 

p_ (Average density of solid) » 1 ,Zh g/cmL C « 3% Bl 

5 - 7.36 g/cmj, C° » 15% Bi 

* y.b6 g/cm , C° * 30% BI 

The density of the solid (tin rich - 0 phase) Is determined using lattice 
spacing for pure Sn as a function of temperature, taken from Handbook of Lattice 
Spaclngs and Structure of Materials (HLSSM). The variation of density with 
temperature agrees exactly with that computed by using the thermal expansion 
coefficient given in Smlthells, p, 9*<1, and there Is excellent agreement with 
the density (26°C) given In Handbook of Chemistry and Physics , 57th ed., CRC 
Press (HCP). The effect of BI on the density of 8 - Sn Is from h'LSSM where the 
lattice spaclngs of 8 - Sn Increase linearly with atom fraction BI. At temp- 
eratures of the solidification range, the same variation is assumed. 

For macrosegregation calculations, the average density of the "cored" solid 
phase Is used; as such, the density of the solid depends upon the average 
composition of the solid. The average composition of the solid Is estimated 
by the Schell equation. 

The density of eutectic solid is computed by using the density of BI at 26°C 
and correcting for temperature wl th the thermal expansion coefficient (these 
data are in HCP) . No data on the effect of Sn on density of Bi are available 
so the density of Bi - rich eutectic solid (98. b% BI) is computed assuming 
Ideal solution; i.e., the molar volume Is given by the "law of mixtures". 

Pens! ty of Liquid 

P L = 0.0277 C L + 6.99 g/cnr 5 

(density of eutectic liquid) = 8.65 g/cm^ 

Density of liquid pure Sn is taken to be the average of data in Thresh and 
Crawley, Smithells and HCP all extrapolated to the solidification temperature 
range. In this range, the data from these sources agree well. For Bi four 
sets of data are in HCP and another in Smithells. Three sets of data which 
showed closer agreement were selected and averaged. However, whether 3 sets 


or all 5 sots of data are considered makes little difference. 

The density of liquid alloys was computed by assuming that the solutions are 
Ideal; l,e», there are no volume changes due to solution and the "law of 
mixtures" applies. The procedure can be justified since the heat of mixing 
for these alloys Is small Indicating almost Ideal behavior. 

Viscosi ty 

y « 1.583 x lO -11 C L + 0.0163 

with y In Poise (g/s~cm) and Is the composition along the 
Hquldus (%Bi). 

If constant values of y are selected: 

y * . 0?. 15 g/s-cm, C q * 3% B1 

» 0.023 g/s-cm, C q ** 15% Bi 

* 0.024 g/s-cm, C » 30% BI 
o 

For pure Sn; Thresh and Crawley show that the viscosity obeys the Andrade 
equation; thus this Is used as a basts to extrapolate the viscosity Into the 
solidification temperature range. Viscosity data (from HCP) for pure Bi are 
reduced to the form of the Andrade equation and also extrapolated to the 
solidification temperature range. Viscosity of alloys is assumed to follow 
the "law of mixtures". 
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